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It is shown that the ratios of the quenched averaged three and four-point correlation functions of 
the local energy density operator to the connected ones in the random-bond Ising model approach 
asymptotically to some universal functions. We derive the explicit expressions of these universal 
functions. Moreover it is shown that the individual logarithmic operators have not any contribution 
to the connected correlation functions of the disordered Ising model. 
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Random systems represent the spatial inhomogenuity 
where scale invariance is only preserved on average but 
not for specific disorder realization. The understanding 
of the role played by quenched impurities of the nature of 
phase transition is one of the significant subjects in sta- 
tistical physics and has attracted a great deal of attention 
[1]. According to the Harris criterion [2], quenched ran- 
domness is a relevant perturbation at the second-order 
critical point for systems of dimension d, when its spe- 
cific heat exponent a, of the pure system is positive. 
Concerning the effect of randomness on the correlation 
functions, it is known that the presence of randomness in- 
duces a logarithmic factor to the correlation functions of 
pure system [3] . Theoretical treatment of the quenched 
disordered systems is a non-trivial task in view of the 
fact that, one has to average the logarithm of the parti- 
tion function over various realization of the disorder in 
the statistical ensemble and therefore find physical re- 
sults. There are two standard methods to perform this 
averaging, the supersymmetry (SUSY) approach, and the 
well-known replica approach. Recently using the replica 
approach it has been shown by Cardy [4], that the log- 
arithmic factor multiplying power law behavior are to 
be expected in the scaling behavior near non-mean field 
critical points. It is shown also that the results are valid 
for systems with short-range interactions and in an ar- 
bitrary number of dimensions. He concludes that in the 
limit of n — > of replicas the theory posses of a set of 
fields which are degenerate (they have the same scaling 
dimensions) and finds a pair of fields which form a Jordan 
cell structure for dilatation operator and derives logarith- 
mic operator in such disordered systems. Cardy proves 
that the quenched disordered theory with Z = 1 can be 
described by logarithmic conformal field theory as well. 
The logarithmic conformal field theories (LCFT) are ex- 
tensions of conventional conformal field theories, which 
have emerged in recent years in a number of interesting 
physical problems of condensed matter physics [5] , string 
theory [6], and nonlinear dynamical systems [7]. The 



LCFT are characterized by the fact that their dilatation 
operator L are not diagonalized and admit a Jordan cell 
structure. The non-trivial mixing between these opera- 
tors leads to logarithmic singularities in their correlation 
functions. It has been shown [8] that the correlator of two 
fields in such field theories, has a logarithmic singularity. 
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In this direction we show that the quenched averaged 
connected correlation functions of local energy density 
field can be written in terms of ordinary scaling operators 
which can be constructed by the difference of energy op- 
erators in two different replicas. We write the connected 
3 and 4-point correlation functions of energy density ex- 
plicitly in terms of such ordinary operators. Further- 
more we prove that the logarithmic operators have no 
contribution in the quenched averaged connected corre- 
lation functions of the local density operator. However, 
these operators play a considerable role on the discon- 
nected ones and produce some logarithmic factors in the 
correlation functions. We calculate the various types of 
quenched averaged 3 and 4-point correlation functions of 
the local energy density and show that the ratios of these 
correlation functions to the connected ones have the spe- 
cific universal asymptotic and write down these universal 
functions explicitly. 

2-Connected Correlation Functions of Random 
Magnets 

We consider a quenched random ferromagnet, for in- 
stance an Ising model, with random-bond disorder. Let 
us describe this disordered system in the continuum limit 
by the following Hamiltonian, 



H = H° 



J(r)E{r)d d r 



(2) 



where H° is the Hamiltonian of the renormalization 
group at fixed point describing the pure Ising model. 
The field J(r) is a quenched random variable coupled to 
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the local energy density E{r). When the coupling J(r) 
is independent of x and not random, the above Hami- 
tonian describes the behavior of the statistical model 
near it's critical point. For simplicity we assume that 
the random variable J(r) is a gaussian variable which 
is characterized by its two moments < J(r) >= and 
< J(r)J{r') >= gd(r — r'). The standard procedure of 
averaging over disorder is to introduce replicas, i.e., n 
identical copies of the same model for which 

n r, n 

Z n = TreM-J2 H «- d d rJ(r)J2E a (r)} (3) 
0=1 J 0=1 

averaging over the disorder gives rise to the following 
effective replical Hamiltonaian. 



n 

o=l J a^b 
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we keep only the non-diagonal terms, since using the op- 
erator algebra of the pure system one can absorb the 
diagonal terms into H®. The replicas are now coupled 
via the disorder. The scaling dimension of coupling g is 
y g = d — 2A E and is relevant at the pure fixed point 
if y g > 0. For small y g it is possible to use standard 
perturbation theory and find the possible random fixed 
point. It is noted by Cardy that the n operator E a are 
degenerate at the pure fixed point and one can decom- 
pose them into irreducible representation of permuta- 
tion group S n . It has been shown that the combination 
E t — X)"=i E a is a singlet ( symmetric under the permu- 
tation of the replica group) and E a = E a — — Y^b=i Eb 
transforms according to an (n — l)-dimcnsional repre- 
sentation of S n . The fields E a satisfy the condition 
S"=i E a = 0. The important observation is that the 
fields E t and E a have the proper scaling dimensions close 
to n -> as A Et = A^ + ±(1 - n)y g + 0(y 2 g ) and 

= A^ + \yg + 0(y 2 ) respectively. It is clear that the 
singlet field E t becomes degenerate with the (n — 1) oper- 
ators E a . This is true to all orders [4]. However they do 
not form the basis of the Jordan cell for the dilatation op- 
erator. To find the logarithmic pair according to [4] we 
define the correlation function of < Et(0)Et(r) >= A\ 
and < E a (0)E a (r) >= B\ and find the following rela- 
tions for Ax and B\. 

A x = n(a - (n - 1)6) = nA{n)r- 2AE ^ 

Bx = (1 - -)(a - b) = (1 - -)B(n)r- 2KE{ ^ (5) 
n n 

where a =< E^OjE^r) > and b =< Ei(0)Ej(r) > 
with i ^ j. the above equations enable us to write the 
quenched averaged connected two-point correlation func- 
tions of energy density operator in terms of a and b in 
the limit of n -> as: < E(0)E(r) > c = a - b which 



is equal to B(0)r~ 2AE and it has a pure scaling behav- 
ior. However, the correlation functions a and b have the 
logarithmic singularities and behave as: 

<Ex{0)Ex(r)> = (A'(0) - B'(0) + S(0) 



5(0)^ In r)r- 
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< Ex(0)E 2 (r) > = (A'(0) - B'(0) - A(0)^lnr) r - 2AE (6) 

where A(0) — B(0). The prime sign in the eq. (6) means 
differentiating with respect to n. This means that in the 
limit n — > the fields E t and E a form a basis of Jordan 
cell, i.e. their two point correlation functions behave as: 

< E t (0)E t {r) >= 0, < E t (0)E a (r) >= ai r- 2AE and 

< E a (0)Eb(r) >= (-2ai lnr + D a . b )r~ 2A , where a\ and 
D a> b are some constants. As noted by Cardy the ratio 
of quenched averaged two-point correlators of the energy 
density operator to the connected one has a universal 
r-dependence as: 



< E(0)E(r) > < £7(0) >< E(r) > 



< E(0)E(r) > c < E(0)E(r) > c 



lnr 



(7) 



To understand the structure of Jordan cell, we note 
that in 2D one can define the operator L as 



Lo 



A E 
1 A B 



(8) 



so that, L E t = A E E t and L E a = A E E a + E t in the 
limit of n — > 0. Using this representation for Lq one 
can show that the field E t with its logarithmic partner 
E a have the standard logarithmic correlation functions 
[9], (see the correlation above the eq.(8) ). We note 
that in 2D we have dealt with two-dimensional confor- 
mal field theory, relying heavily on the underlying Vi- 
rasoro algebra. For an extention to D dimensions one 
has to modify the representation of the Virasoro alge- 
bra to higher dimensions [10]. We consider a doublet of 

fields (Jordan cell) $ = ^ ^* ^ and note that under D- 

dimcnsional conformal transformation x — > x', we have, 
$(x) — ► $'(x') = G T $(x) where T is a two dimensional 
matrix which has Jordan form and G = ||^-|| is the 
Jacobian. For our particular case T has the following 
Jordan form: 



T= [ x ~ U 



2Ae 
D 



(9) 



and one can show that the two fields E t and E a , trans- 
form as: 



E t (x') = G-^E t ( X ) 

K(x') = G ~ (]n(G)E t (x) + E a {x)) 



(10) 
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This expresses that the top-field E t always transforms as 
an ordinary scaling operator. It can be verified that the 
correlation functions of fields Et and E a have the stan- 
dard D- dimensional logarithmic conformal field theory 
structure [9,10]. Using the above results, it is evident 
that the dimension of field-difference E a — E b with a ^ b 
is and it transforms as an ordinary operator under 
the scaling transformation. The intresting observation 
is that the connected averaged correlation functions de- 
pends on the difference fields E a — E b only and therefore 
they behave as the ordinary correlation functions. For in- 
stance in the following we write the connected quenched 
averaged 2,3 and 4-point functions of local energy density 
in terms of the field-difference operators explicitly, 

1 



< E(1)E(2) >c = -<( Ea - E b ) (1) (E a - E b ) (2) > (11) 



< E(l)E(2)E(3) > c = < (E a - 

{E a -E c ) {2) {E a -E b ) m > (12) 



< E(1)E(2)E(3)E(A) > c = 



< {E a - E b )^(E a - E c )(2){E a ~ E d )(3){E a - £*)(4) > 



— 2 < (E a - E b )(i)(E c - E d )( 2 ){E c - Ed)(3)(E a - E b )^ > 
< {E a - Eb)(i)(E c - E d )( 2 ){E a - Eh)(3)(E c - E d )^ > 

-\ < (E a - E b ) (1) (E a - E b ) [2) (E c - E d ) (3) (E c - E d ) (i) > 

where the last equation has only 15 independent terms. 

To confirm this prediction also one can directly show 
that the quenched averaged connected correlation func- 
tions have a pure scaling behaviour which is determined 
by ordinary scaling operators and the logarithmic opera- 
tors E a do not change its behavior. This can be verified 
directly for the quenched averaged connected 3-point cor- 
relation function of energy density. 

We are interested in deriving exactly the various 3- 
point quenched averaged func- 

tions as < E(1)E(2)E(3) >, < E(l)E(2) >< E{3) > 
and <E(1)>< E(2) >< E(3) >, which can be writ- 
ten in terms of the replica correlation functions < 
Ei(l)E! (2)£i (3) >= a < E 1 (1)E 1 (2)E 2 (3) >= b and 
< Ei (l)E2 (2)E% (3) >= c, respectively. One can derive 
the correlation functions a, b and c by means of 3-point 
functions of E t and E a as follows: 



< Et(l)E t (2)E t (3) >= na + 3n(n - 1)6 
n(n — l)(n — 2)c = nA\ 



< E a (l)E a (2)E t (3) >= n ia + (n?(n - 1) 



_4n 2 + -^(n-l) 2 + -|(n-l)(n-2))6 
+ (~ ni (n - l)(n - 2) + -L(n - 2) 2 (n - l))c 
= (1 - 

n 



(14) 



and finally, 



< E a (l)E a (2)E a (3) >= (nf - ^^)a 



n-1. 



" ~ 1)(n ~ 2) + h niin ~ l))h 

+(ini(n - l)(n - 2) - \{n - l)(n - 2)(n - 3))c 

B (l-I)(l-|)Ci (15) 

where n\ = (1 — ^) and A 1; i?! and C\ are pure scaling 
functions of variables r^j . To derive the above equations 
we use the replica symmetry and symmetry of the various 
types of 3-point correlation functions under interchang- 
ing of positions. We note that replica symmetry leads 
to have < E a {l)E t (2)E t (3) > = and therefore, dose 
not give any new relationship between a, b and c. Using 
the above equations, it can be found that the correlation 
functions a, b and c are as follows: 



a = 



3nBi - 3nCi + n 2 Ci + Ai - 3Bi + 2Ci 



nB 1 - nCi +A 1 - 3B X + 2C X 



A 1 - 3B X + 2Ci 



(16) 



Using the above equations we can show that the con- 
nected quenched averaged 3-point function behaves as: 



< E(1)E(2)E(3) > = 2c + a-3b 



Ci 



(17) 



(13) 



which is a scaling function and confirms the observa- 
tion that the logarithmic operators (individually) have no 
role in the connected quenched averaged correlation func- 
tions. In addition one can derive the correlation functions 
< Ei(l)Ej(2)E k {3) > for given i, j and k in the limit of 
n — > and show that they have the following form: 

< £,(1)^(2)^(3) > = [a ijk - PijkDx 

+ 7 y - fe (4^2-^ 2 )]/(l,2,3) (18) 
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where /(1,2,3) = (r 12 r 13 r 23 ) 2Ab , D 1 = ln(r 12 r 13 r 2 3) 
and D 2 — In r 23 In r\ 3 + In r\ 3 In r\ 2 + In In r\ 2 . It can 
also be shown that the ratio of various symmetrized 3- 
point functions to the connected one behaves asymptot- 
ically as a universal function 



~(4£> 2 -£>?) 



(19) 



We generalize the above calculations to derive the var- 
ious type of 4-point correlation functions and show that 
the ratio of the various disconnected to the connected 
one have the following universal asymptotic: 



60 2 - 30 3 - 120 4 - I8O5] 



(20) 



where Oi = ln(r 12 r 13 ri4r 23 r 2i r 3i ), 2 = (lnr^ lnrjy + 
■ • •) with i^j^k^l,0 3 ~ (lnry lnr^, + ■ ■ •) with i ^ 

j ^ k, Oi — (bxrij lnr fc; lnry H ) with i j k I, 

and finally O5 = (In In nk In ru + • • ■) with i j 
k^l 

In summary, in this paper we have studied the correla- 
tion functions of disordered random magnets [11] and ob- 
tain the various types of 3 and 4-point quenched averaged 
correlation functions. One can check directly that these 
different types of the 3 and 4-point correlation functions 
have the general property of a logarithmic conformal field 
theory that the logarithmic partner can be regarded as 
the formal derivative of the ordinary fields (top field) with 
respect to their conformal weight [9]. In this case, one can 
consider the E a fields as the derivative of E t with respect 
to n . We emphasise that the derivative with respect to 
scaling weight can be written in terms of the derivative 
with respect to n. These properties enable us to cal- 
culate any N-point correlation function containing the 
logarithmic field E a in terms of the correlation functions 
of the top-fields. The general expression of the correla- 
tion functions of the LCFT's are given in ref. [9] and here 
we determine the unknown constants in the logarithmic 
correlation functions in terms of details of the random- 
bond Ising model. It is noted that the formal derivations 
with respect to scaling dimensions can not predict the 
unknown constants in the quenched averaged correlation 
functions of the local energy density operators. The con- 
stant depends on the detail of the statistical model. We 
have shown that the individual logarithmic operators E a 
do not have any contribution to the quenched averaged 
connected correlation functions of the energy density. We 
also obtain that the connected correlation functions can 
be written in terms of the difference fields which trans- 
form as an ordinary scaling operator. However they will 
play a crucial role to the disconnected averaged correla- 
tion functions. Also we find that the ratio of the vari- 
ous types of 3 and 4-point quenched averaged correlation 
functions to the connected ones have a universal asymp- 
totic behavior and give their explicit form. These predic- 
tions can also be investigated numerically. Our analysis 



are valid in all dimensions as long as the dimension is 
below the upper critical dimensions. To derive the above 
results we have used the replica symmetry. Any attempt 
towards the breaking of this symmetry will change com- 
pletely the above picture and produces more than one 
logarithmic fields in the block and produces higher order 
logarithmic singularities [9]. 

These results can be easily generalized to other prob- 
lem such as polymer statistics, percolation and random 
phase sine-Gordon model etc. 
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